Optimal fast single pulse readout of qubits 
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The computer simulations of the process of single pulse readout from the flux-biased phase qubit 
is performed in the frame of one-dimensional Schroedinger equation. It has been demonstrated that 
the readout error can be minimized by choosing the optimal pulse duration and the depth of a 
potential well, leading to the fidelity of 0.94 for 2ns and 0.965 for 12ns sinusoidal pulses. 
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In a past decade a serious progress has been achieved in 
development and creation of various circuits for quantum 
computation [H,0]. Different sources of decoherence is, 
however, the main factor, limiting the practical utiliza- 
tion of complex networks of quantum bits [^|-[lo|. Re- 
cently, it has been demonstrated that, the coherent Rabi 
oscillations remain nearly unaffected by thermal fluctua- 
tions up to temperatures of IK [ll| (i.e., until the energy 
of thermal fluctuations kT becomes comparable with en- 
ergy level spacing hu> of the qubit), so without degrad- 
ing already achieved coherence times, phase qubits can 
be operated at temperatures much higher than those re- 
ported till now. This may signal, that relatively large 
readout errors of practical devices can be attributed not 
to quantum and thermal fluctuations, but to unoptimal 
readout of the qubits. To speed up the readout, the fast 
single pulse readout (FSPR) technique has been realized 
and tested 0-0 ■ An example of a shallow potential well 
with the two energy levels |0) and |1) is presented in the 
inset of Fig. 1. The basic idea of the FSPR is the ap- 
plication of a readout pulse in such a way, that during 
the pulse action the system will tunnel from the state 
|1) through the barrier with a probability close to unity, 
while from the state |0) it will not tunnel, again, with 
the probability close to unity. The effect of the shape 
and duration of the pulse on different error probabilities 
has been studied in Ref.s , [lj| • In particular, the error 
to excite higher qubit states due to nonadiabaticity of 
the pulse was analyzed [loj . As it has been understood, 
j [lot > the main source of error during qubit readout is 
due to incomplete discrimination between the two quan- 
tum states |0) and |1). It has been claimed that the 
quality of tunneling discrimination depends significantly 
on the measurement pulse amplitude and this error de- 
creases for longer measurement pulses and with increase 
of Ti/Tq ratio (where Tq and T% are tunneling rates from 
the states |0) and |1), respectively). However, the ques- 
tion about possible compromise between the speed and 
the fidelity of the readout has not been considered in the 
literature. The importance of this question is not only 
due to the reason that one wants both high fidelity and 
high speed readout, but also due to the fact that any fluc- 
tuations will accumulate during longer readout, leading 
to entangling of states. 
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FIG. 1: The profile of a bistable potential. Dashed curve - 
the original potential, solid curve - the potential with enlarged 
deep well to simulate the effect of damping. The inset: the 
enlargement of a shallow potential well. 



Recently, in classical systems subjected to noise and 
pulsed or periodic driving, such as Josephson junctions 
12] and magnetic nanoparticles [l3j |. it has been demon- 
strated that at a fixed value of the driving amplitude 
there exists the optimal pulse duration, which minimizes 
the noise-induced errors. The aim of the present paper 
is to study the readout error N, which is a sum of the 
two probabilities, Pro not to tunnel during the pulse ac- 
tion from the state |1), and Poi to tunnel from the state 
|0) (i.e. N = Pro + P i, while the fidelity F = 1- N). 
It is intriguing to understand, would the effect, discov- 
ered for classical systems with noise 12j, 13j, be real- 
ized in a quantum system, e.g. in a such, as an example 
of a qubit, described in , [3] , [n| . The investigation is 
performed via computer simulation of the Schroedinger 
equation and is focused on the dependence of the readout 
error N versus the pulse amplitude, duration and shape, 
as well as the depth of a shallow potential well. 



Let us consider an example of a flux-biased phase qubit 
0,0! which is described by the following poten- 
tial V(x,t) = Ej {{x - <p(t)) 2 /2£ - cos a;}, see Fig. 1, 
dashed curve. Here Ej — Ich/2e is the Josephson en- 
ergy, x is the Josephson phase, e is the electron charge, 



and h is the Planck constant. We take for the qubit 
the same parameters as in [6], 10j: the critical current 
Ic = 1.7/iA, the inductance of the ring L — 0.72nH 
and the capacitance C = 700fF, thus I = 2eIcL/Ti = 
3.71, 2e 2 /hC = 0.6933 x 10 9 Hz, Ej/h = I c /2e = 
5.31 x 10 12 Hz, so it is convenient to introduce the "in- 
verse capacitance" D = 2e 2 /hC x 10~ 9 and express the 
time in nanoseconds. The dimensionless external mag- 
netic flux <f(t) consists of the two components, the dc 
component a , adjusting which the depth of the shal- 
low well can be changed, and the driving readout pulse: 
ip(t) = 2n(ao + Af(t)), where A is the pulse amplitude, 
and f(t) is one of sin(7rf/i p ), sin 2 (nt / 't p ) , sin 4 (7rf/i p ), the 
trapezoid function, which linearly grows and drops for 
t > tp/4 and t < 3t p /4, and the sine-trapezoid func- 
tion, which differs from the previous one by the sinu- 
soidal walls, see the inset of Fig. [51 We note that t p is 
defined as the full width of the pulse at zero level, not 
as full width at half maximum. The shift of a potential 
barrier ao = 0.81 is chosen such to allow the six levels 
to be present in a shallow potential well, see [loj | . The 
pulse with the amplitude A rs 0.035 leads to lowering the 
potential barrier such that only two levels will remain. 

The Schroedinger equation for the wave function 
^>(x,t) has the following form: 



,dV(x,t) 2e 2 <9 2 *(x,i) V(x,t) 



dt hC dx 2 h 

It is assumed, that the boundary conditions 

V(c,t) = $(d,t) = 0, 



*(*,*). (1) 



(2) 



at the points c and d are taken far away from the shal- 
low potential well, and do not affect the tunneling pro- 
cess. To prevent from the repopulation error 10], which 
arises due to the absence of damping in our model, let 
us introduce an effective damping in the following man- 
ner. Since we are interested in the process of tunneling 
from the shallow potential well only, we assume that at 
the bottom of the deep potential well the potential does 
not grow up, but continues to the right really far away, 
see Fig. 1, solid curve. In particular, we have taken 
c = —3, d = 797, while the left minimum is located 
at xi « 1.4, and the right one at X2 ~ 6. Numerical 
solution of the Schroedinger equation (fTJ) with boundary 
conditions (J2J) has been performed on the basis of implicit 
finite-difference Crank-Nicholson scheme. Typical values 
of discretization steps are Ax = 0.01, At = 10~ 4 ns. 

In Fig. 2 the evolution of the probabilities P(t) not 
to tunnel from the state |1) and Q(t) to tunnel from the 
state |0) are presented. In both cases we choose the initial 
condition (the wave function) to correspond to either |1) 
or |0) stationary state and calculate the survival proba- 
bilities Pi (i) and Pq (t) in the shallow potential well. The 
probability evolutions P(t) = P x (t) and Q(t) = 1 - P Q (t) 
are presented in Fig. 2 for the case of sinusoidal pulse 
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FIG. 2: The evolution of probabilities P(t) and Q(t). In- 
set: the considered pulses, from top to bottom, sine-trapezoid 
pulse, trapezoidal pulse, sm(irt/t p ), sin 2 (-Kt/t p ), sm 4 (-irt/t p ). 



with driving amplitude A — 0.034 and pulse duration 
t p = 8ns. One can see that at the end of the pulse both 
probabilities are significantly smaller than unity, and our 
task is to find the optimal parameters of the pulse to min- 
imize the readout error N = P w + P i = P{tp) + Q(tp)- 
In difference with [l(| we do not separate the nonadia- 
batic error (leading to lifting to higher states), and the 
error of direct tunneling from the state |0), because both 
these errors contribute into the error Pox , which we calcu- 
late. Since for short pulse durations the probability Pio 
must be large, and for long pulses the probability Poi 
increases significantly, there must be an optimal pulse 
duration, leading to the minimal N. 
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FIG. 3: The readout error iV versus the pulse amplitude A 
for the pulse sin(nt/t p ). 



The readout error N versus the pulse amplitude A for 
the sinusoidal pulse is presented in Fig. 3 for different 
values of the pulse width t p . It is seen that N is very 
sensitive to the pulse amplitude, especially for large pulse 



2 



durations, and that the minimal value of N decreases 
with increase of t p . Similar figure for the only one pulse 
width has been presented in Ref. [8J]. However, from 
Fig. [3] one can see that for the pulse durations t p > 8ns 
the minimal value of N is almost the same, while below 
t p < 2ns N increases significantly. 
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FIG. 4: The readout error N versus the pulse width t p 
for the pulse sin(7ri/i p ). The dashed curves correspond to 
sin 2 (nt/t p ), sin 4 (7rf/t p ) pulses with A — 0.036 (from left to 
right). 

The readout error TV versus the pulse width t p for the 
sinusoidal pulse is presented in Fig. [4] for different val- 
ues of the pulse amplitude A. One can see that N has 
strongly pronounced minimum, and with increase of A 
this minimum shifts to smaller tunneling times, corre- 
spondingly speeding up the readout. With increase of A 
the value of the error N at the minimum increases, so one 
should find the compromise between the error and the 
speed. However, in the range of small amplitudes N in- 
creases insignificantly, compare the curves for A = 0.03, 
A = 0.032 and A — 0.034, and for the considered pa- 
rameters the amplitude A — 0.034 can be chosen as a 
compromise value, leading to the fast readout with the 
high fidelity F — 0.94. The results, presented in Fig.s 
[3] and |4j on one hand, demonstrate similar dependence 
of the readout error vs pulse width as for classical sys- 
tems [Hj]) [HI- On the other hand, the dependence on 
the driving amplitude also has a minimum, which out- 
lines the quantum nature of the described phenomena, 
while for a classical system the error monotonously de- 
creases with increase of the pulse amplitude. Another 
important deviation from a classical system is the impos- 
sibility to use the rectangular readout pulses (which in a 
classical case leads to the minimal noise-induced errors). 
In the qubit the rectangular pulse leads to nonadiabatic- 
ity of the tunneling event that results to lifting to higher 
eigenstates and considerable increase of the probability 
to tunnel from |0) state and thus to much larger values 



of N than for all other considered pulse shapes. 
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FIG. 5: The readout error N versus the pulse width t p 
for the trapezoid pulse with A = 0.034 (thick solid curve) 
and sine-trapezoid pulses (solid curves with triangles, dia- 
monds and circles for amplitudes A = 0.034; 0.032; 0.03). 
The curves for the pulse shape sin(nt/t p ) from Fig. [4] with 
A = 0.038; 0.036; 0.034; 0.032 are given by dashed curves for 
comparison. Solid curve with crosses - long sine-trapezoid 
pulse for A = 0.03. 

The readout error for different pulse shapes sin 2 (7ri/t p ) 
and sin 4 (7rf/t p ), A — 0.036 is presented in Fig. 2] by 
dashed curves. As one can see, the minimal value of the 
error is pretty much the same, but the optimal readout 
time is significantly smaller for sin(7rf/i p ) just due to the 
fact, that the top of the pulse in this case is more flat. In 
accordance with these results the trapezoid pulse should 
lead to smaller N. This is indeed so, as one may see 
from Fig. [5j thick solid curve. However, the dependence 
of N vs t p demonstrates oscillations, and it is not easy to 
properly choose the optimal pulse width. Better situation 
can be achieved, if one substitutes the linear walls in the 
trapezoid by sin(27ri/i p ), see solid curves with triangles, 
diamonds and circles for amplitudes A — 0.034, A — 
0.032 and A = 0.03, respectively. In comparison with 
pure sinusoidal pulse, sin(7rt/t p ), one can gain about 15- 
20 % in the readout error, and, simultaneously, up to 
50 % in the readout speed, compare with the dashed 
curves, taken from Fig. H for A = 0.038, A = 0.036, 
A = 0.034, A = 0.032, respectively. Further increase 
of the readout speed can be achieved, if one increases 
the duration of the flat part of the sine-trapezoid pulse 
from tp/8 to 7t p /8, see solid curve with crosses for A = 
0.03. Transferring to the limit of a rectangular pulse, one 
should stop somewhere, since further increase of the flat 
part of the pulse will give little increase in speed, but will 
obviously lead to the increase of nonadiabatic error [l(| • 

Finally, let us briefly demonstrate, how adjusting the 
depth of a shallow potential well, the readout error can be 
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FIG. 6: The readout error N versus the pulse width t p for 
the pulse sin(7r£/t p ) with different amplitudes and inverse ca- 
pacitances D. Inset: the readout error N(D) for the pulse 
duration t p = 12ns (circles) and t p = 2ns (diamonds). 



pulse, whose generation does not require complex pulse 
shaping hardware (l4| . For 2ns sine-trapezoid pulse the 
fidelity can be increased to 0.947. We, therefore, believe 
that the proper engineering of the qubit, together with 
the optimal adjustment of the readout pulse duration, 
will finally lead to creation of high fidelity and high-speed 
readout qubits. 

In conclusion, we have performed the computer simula- 
tions of the process of pulsed readout from the flux-biased 
phase qubit within the model of one-dimensional time- 
dependent Schroedinger equation. It has been demon- 
strated that by choosing the optimal pulse duration the 
readout error can be minimized. Further decrease of the 
readout error can be achieved by variation of the depth 
of the shallow potential well. 

The authors wish to thank M.A. Silaev for discussions 
and useful comments. The work has been supported by 
RFBR (projects 09-02-00491 and 08-02-97033). 



decreased further. If this well is deep enough, there are 
many levels inside, and while both Tq and T\ are small, 
their values are close to each other, which complicates 
the discrimination between the states |0) and |1). If the 
well is too shallow, the tunneling from the state |1) may 
occur even without the driving pulse, and the discrimina- 
tion between the states, again, will be pure. Therefore, 
there must be the optimal depth of the well, leading to 
the minimal readout error. The depth of the well can 
be adjusted either by variation of the constant magnetic 
field component ao, or by variation of the capacitance C 
of the Josephson junction. The latter case is illustrated 
in Fig. [6l where the readout error N is presented for dif- 
ferent amplitudes and values of the inverse capacitance 
D = 2e 2 /hC x 10~ 9 in such a way that the minimum 
of N remains in approximately the same region of pulse 
duration t v rs 12ns. Fixing t p , from the main part of 
Fig. [5] one can extract the values of N and plot them 
as a function of D, as it is done in the inset of Fig. [6l 
the curve with circles. The readout error N(D) demon- 
strates pronounced minimum and choosing the optimal 
value of D rj 1.4Hz one can decrease the value of N down 
to 0.0355, leading to the fidelity F = 0.965 (note that, 
N = 0.0523 for D = 0.6933Hz, A = 0.033). The same 
procedure can be performed for t p = 2ns (see the curve 
with diamonds in the inset): while for D = 0. 6933Hz the 
fidelity was about 0.9, at the minimum, which is shifted 
to D w 1.8Hz, F ks 0.94. The location of the minimum 
of N(D) corresponds to something between three or two 
energy levels inside the shallow potential well. We note, 
that rather fast readout of 2ns duration with acceptably 
high fidelity of 94% is achieved with a simple sinusoidal 
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